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Scaling regimes of a semi-flexible polymer in a rectangular channel 
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We derive scaling relations for the extension statistics and the confinement free energy for a semi- 
flexible polymer confined to a channel with a rectangular cross-section. Our motivation are recent 
numerical results [Gupta et al., JCP 140 (2014) 214901] indicating that extensional fluctuations 
are quite different in rectangular channels compared to square channels. Our results are of direct 
relevance for interpreting current experiments on DNA molecules confined to nano-channels, as 
many experiments are performed for rectangular channels with large aspect ratios while theoretical 
and simulation results are usually obtained for square channels. 
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Conformations of semi-flexible polymers confined to 
narrow channels are presently intensively investigated, 
since it has emerged that nano-channel confined DNA 
molecules offer possibilities both for biological applica¬ 
tions and as a model system for experimental polymer 
physics [1-10]. These studies, and the theoretical and 
simulations results which have accompanied them, have 
shown that semi-flexible polymers such as DNA exhibit 
a much richer spectrum of behaviours under confinement 
than do flexible polymers [11-16]. 

Since it is convenient to fabricate channels with a fixed 
height but varying widths, many experimental studies [1— 
5, 8, 17-19] of confined DNA are performed in rectangu¬ 
lar channels (width D\ y, height Du) with aspect ratios 
far from unity, D\ y Du- However, most simulation 
and theoretical studies [5, 12-16, 20-22] are restricted to 
channels with square cross-sections, D^j = Du- A com¬ 
mon procedure is to analyse the experimental results in 
terms of the ‘effective channel size’ \/DyjDu, simply dis¬ 
regarding the influence of the aspect ratio. 

Recently Gupta et al. [8] have shown that the variance 
of the extension of the DNA molecule does depend on 
the aspect ratio, but a theoretical explanation for this 
intriguing result is lacking. This motivated us to analyse 
how the aspect ratio influences the extension statistics 
of confined semi-flexible polymers. At first sight this ap¬ 
pears to be a difficult problem, because it is governed 
by a large number of length scales: the contour length 
L of the polymer, its effective width w, its Kuhn length 
Ik [23], its global persistence length g [11], the typical 
contour-length separation between intra-chain collisions 
fee, the typical contour-length separation between colli¬ 
sions with the floor and the ceiling of the channel l c h, 
the typical contour-length separation between collisions 
with the vertical walls of the channel l cw , and the chan¬ 
nel width D\y and height Du- There are many different 
confinement regimes corresponding to different orderings 
of these length scales, potentially resulting in a very com¬ 
plicated phase diagram (Fig. 1). Little is known about 
this phase diagram for rectangular channels, except in 
the limit of very strong confinement [11, 24, 25]. To in¬ 
terpret a given experiment it is necessary to determine 
which regime in the phase diagram it corresponds to, 


and what the resulting scaling relations for the extension 
statistics are. 

In this paper we summarise the results of our analysis, 
based on the mean-field theory for the extension of an 
unconfined semi-flexible polymer [26]. It is well known 
how this theory must be adapted to describe the exten¬ 
sion of semi-flexible polymers confined to square chan¬ 
nels with Z?w = Du = D [27]: for a wide channel and a 
long polymer, the polymer globule is divided into a series 
of smaller spherical blobs of size D. One assumes that 
mean-field scaling holds for each blob and concludes that 
the extension of the confined polymer scales as D ~ 2/, ‘ 3 . 
We generalise this analysis to rectangular channels. To 
this end it is necessary to consider a hierarchy of blobs 
(inset of Fig. 1). For the special case of a flexible chain 
this approach was used by Turban to compute its exten¬ 
sion in a rectangular channel [28]. We emphasise that 
this is a much simpler problem since a flexible chain ex¬ 
hibits only a single confinement regime, as opposed to 
the semi-flexible polymer, as pointed out above. 

An important result of our analysis is that the scal¬ 
ings can be simply summarised, despite the fact that the 
phase diagram Fig. 1 exhibits many different regimes. 
First, we find that the average extension is approximately 
independent of channel aspect ratio, provided that at 
least one of the channel dimensions is significantly larger 
than the Kuhn length of the polymer. This is an impor¬ 
tant finding because it implies that it is approximately 
correct to analyse the average extension of nano-confined 
DNA molecules in terms of the effective channel size 
y/D\yDu- Second, we find that the extension variance 
depends strongly on both channel dimensions separately, 
it would be incorrect to analyse it in terms of the effective 
channel size. We find that the variance increases rapidly 
with aspect ratio if far from unity. Our theoretical results 
for the variance explain the findings of Ref. [8], they also 
make it clear that square channels are much preferred for 
applications where extensional fluctuations are required 
to be as small as possible. Third, we compute the free 
energy of confinement, it is approximately determined by 
the smallest confining dimension. These scaling predic¬ 
tions were derived under the mean-field approximation of 
Flory [26], but in certain parameter regimes (regimes Ha- 
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FIG. 1. Phase diagram of different scaling regimes for the ex¬ 
tension statistics of a semi-flexible polymer confined to a chan¬ 
nel with a rectangular cross-section. The scaling in regime I is 
given by Eqs. (7,8). The scaling in regime II is given by Eqs. 
(16,17). The vertical grey lines distinguish different scaling 
regimes of the confinement free energy, see text. Results for 
regime III have been derived elsewhere [11, 16, 24, 25] and are 
discussed in the SM. Inset: An illustration of the hierarchy 
of blobs analysed in regime I. In regime la, the smaller blobs 
are spherical blobs of size Rbiob = Du- In regime Ib-c they 
are cylindrical, of height Du and width Rbiob > Du- 

b in Fig. 1), the results are supported by an asymptoti¬ 
cally exact theory that was developed for square channels 
[22], but can be generalised to rectangular channels. 

We have summarised these results in a table in the 
supplemental material (SM) [29], also including results 
for very strong confinement that were derived by other 
authors. For other observables (e.g. dynamics, statistics 
of circular polymers, probability of knot formation), the 
scaling properties remain to be determined, but the dis¬ 
tinctions between the regimes we derive here must also 
apply to them. 

Method. We assume that the polymer is described by 
the self-avoiding worm-like chain model [23] with Kuhn 
length <?k and effective width w < £k- The centre-line 
of the polymer is confined to a rectangular channel of 
height Uh and width > -Dh- We analyse the confor¬ 
mational fluctuations of the confined polymer in terms 
of a hierarchy of blobs assuming that self-avoiding ran¬ 
dom walks obeying Flory’s mean-field theory [26] in two 
and three spatial dimensions. To summarise this theory 
assume that the polymer consists of N independent seg¬ 
ments of length £ and excluded volume v. In d = 2, 3 
dimensions the mean-field result for the average of the 
extension R of the polymer globule is [26]: 

R « ( N 3 £ 2 v ) *** . (1) 

Here the symbol ~ indicates a scaling relation, ignor¬ 
ing factors of order unity. The variance of the extension 
is determined by self similarity. The self-avoiding poly¬ 
mer in two or three dimensions has only one macroscopic 
length scale (R) [23]. Thus, the standard deviation of 
the extension must scale as the average extension [30]: 

c trxsR » (N 3 £ 2 v) . (2) 


As we show below, the analysis of the extension statis¬ 
tics of the confined polymer must proceed by different 
steps, depending on the relation between the contour 
length scales Ik- g , Z c c, ich> lew, defined in the introduc¬ 
tion. The global persistence length g only differs appre¬ 
ciably from the Kuhn length Ik at very strong confine¬ 
ment (regime III in Fig. 1). Theories for this regime have 
been derived elsewhere [11, 16, 24, 25], and are briefly 
discussed in the SM [29]. 

Extension statistics. Consider the separation of scales 

£k < l c c < Ich < lew < L . (3) 

How do the average extension and its standard deviation 
depend upon D w and .Dr? Since l cc <C Z c h, the poly¬ 
mer exhibits three-dimensional Flory scaling before its 
first collision with the channel walls. The first collision 
with the ceiling or the floor must occur when a section 
of contour length Z c h has formed a spherical Flory blob 
of diameter Z?r. Applying Eq. (1) with N = l c h/£u and 
v « I'k'uj [31] yields 

Du « (l 3 ^v/£k) 5 Z c h ~ Ph/(4™)]^ 3 • (4) 

The resulting blobs perform a two-dimensional self¬ 
avoiding walk until they have formed a circular ‘su¬ 
perblob’ of diameter % (we follow the terminology of 
Ref. [28]). This two-dimensional random walk is illus¬ 
trated in the inset of Fig. 1. The number of small spher¬ 
ical blobs constituting one superblob can be estimated 
from Eq. (1) with d = 2, and £ = Dr since each small 
spherical blob constitutes an independent segment of the 
walk. We assume that the small blobs do not overlap, 
thus the two-dimensional excluded volume of the ran¬ 
dom walk of small blobs approximately equals the area 
of a circle of diameter Z)r , in other words v ~ D ^ in this 
context. Eq. (1) results in 

£>w ~ (-^blobs^H) 4 O -^blobs ~ (D W /Dh) 4 D . (5) 

The contour length stored within a superblob equals Z cw , 

lew ~ -^blobs^ch ~ Pw^h/(4»)]^ 3 ' (6) 

There are L/Z cw superblobs that line up along the chan¬ 
nel [Fig. 1 (inset)]. Each superblob has average diameter 
Z?w with fluctuations of the same order. The average 
extension and its fluctuations are therefore given by: 

R ss (L/Z CW )D W ss L [£ k w/(DhDw)] 1 / 3 , (7) 
4 * (T/Z cw )^w » L (^wD^/Dk) 1/3 • (8) 

Eq. (7) was recently derived in Ref. [32], but without 
specifying under which conditions the derivation is valid. 
We now answer this question. The inequality Z cc -C Z c h 
requires that ideal scaling within a blob of size Hr results 
in a large number of intra-chain collisions within the blob 

(.Dr/£ k ) VT>r » D u w/£ 2 k » 1 ^ D h > 4 /w. (9) 
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This corresponds to regime la in Fig. 1. 

The inequality l cw l c h is satisfied if D>w ^ -Dh- 

However, the scaling relations (7,8) reproduce exactly 
the well-known relations for square channels as the limit 
D>w -Oh is taken, commonly referred to as ‘de Gennes 
scaling’ [27] . Thus de Gennes scaling for square chan¬ 
nels is simply a special case of the more general scaling 
relations derived above for rectangular channels. Yet at¬ 
tempting to generalise from square to rectangular chan¬ 
nels by simply replacing D by the geometrical average 
(.Dh-Dw) 1 / 2 gives the wrong prediction for the variance 
of the extension, as Eq. (8) shows. 

Now consider a different ordering of length scales: 

£k < kh < Ice « lew « L . (10) 

This corresponds to Dh <C £\ /w while Dw remains 
large, regime lb in Fig. 1. Are the scaling relations (7) 
and (8) modified in this regime? Since l c h <C l c o the 
polymer must exhibit ideal scaling even after the first 
collision with the ceiling or the floor of the channel. This 
scaling persists approximately until the first intra-chain 
collision, which occurs after a contour length l cc . This 
length scale is estimated by assuming that such a section 
of the polymer will form a cylindrical blob with height 
Z?h and diameter i?biob ~ V^cc^Kj and by setting the 
expected number of collisions within the blob to unity: 

(W^k) 2 ^k w _ l ccW _ . _ 47>h , . 

e >2 71 — D n ~ t t cc — . (it) 

-^blob-^H ZkDh W 

These blobs perform a two-dimensional self-avoiding walk 
until the first collision with the side walls [Fig. l(inset)]. 
Flory scaling for the N = l cw /l cc blobs of size i « i?biob ~ 
V^cc/k and excluded area v « 4iob ~ Icc£k yields 

l cw « [D^D k /(£kw)} 1/3 • (12) 


are identical for these regimes. But note that the scal¬ 
ing analysis is formulated in terms of blobs that obey 
Flory’s mean-field scaling. As long as mean-field theory 
is used throughout, the final prediction for the extension 
must be independent of the way in which the blobs are 
arranged. Since in mean-field theory the repulsive effect 
of self-avoidance is directly determined by the number of 
monomers within the volume spanned by the polymer, 
this also explains why the scaling of the extension is a 
function of the cross-section only, independent of the as¬ 
pect ratio of the channel. 

What about the variance of the extension? Combining 
mean-field theory with the universality of self-avoiding 
random walks [30] shows that the contour length con¬ 
tained in one superblob [Fig. l(inset)] is identical for 
these regimes, and that each blob experiences size fluc¬ 
tuations of order Dw. The variance of the extension 
is given by summing the variances of each blob, yield¬ 
ing a% « fV superb iob s -Dw ~ L(£ K wD%y/D K ) 1/3 . Thus 
the variance increases as the aspect ratio increases. This 
demonstrates that rectangular channels are in fact quite 
different from square ones in this regime. 

If the excluded volume of the polymer is so small 
that the polymer experiences multiple collisions with side 
walls, floor, and ceiling between each intra-chain colli¬ 
sion, then we obtain different scaling relations. Consider 
the following ordering of length scales 

£k < lc h < lew (14) 

It corresponds to Ik <C Dh < Dw <C (Du^k/ w ) 1 ^ 2 
(regime Ha in Fig. 1). Under these conditions the poly¬ 
mer obeys ideal scaling until a blob forms that fills the 
channel cross-section, and is further elongated along the 
channel direction, until it reaches an extension -Rbiob ~ 
ylccfe- The length scale l cc is estimated in a similar way 
as Eq. (11) is obtained: 


This is the same as in regime la, Eq. (6). From this point, 
the derivation of the extension statistics follows that of 
regime la. We infer that the scalings Eqs. (7,8) hold also 
in regime lb. 

Further decreasing Du below Ik one enters a differ¬ 
ent regime (labelled lc in Fig. 1). It corresponds to this 
ordering of length scales: 

lc h < £k < lcc < lew < L . (13) 

In this case the polymer runs almost completely parallel 
to the floor of the channel but can otherwise rotate freely. 
While the expression for the excluded volume of a Kuhn 
length segment must change because of confinement, the 
scaling v « £^w still holds [31]. The same analysis as 
in regime lb can be carried out here, leading to identical 
scaling predictions for the extension statistics, Eqs. (7,8). 
The prefactors will differ between the regimes, however. 

Since the steps in the above derivation are different 
between regime la and regimes Ib-c, it is at first glance 
surprising that the scaling predictions for the extension 


(Icc/£k) 2 £ 2 k w _ ~ / Ik-Ph-Pw \ 1/3 

-Rblob-D H -D W ~ cc ~ V w 2 ) ' 


(15) 


The polymer arranges itself into L/l cc blobs of size i?biob ; 

R « (L/l cc )R Uoh « L [£ k w/ (D H -D w )] 17 f (16) 

4 « (P/U^biob » L£ k . (17) 


We see that the scaling of the average extension agrees 
with (7). The reason is that both equations were derived 
assuming mean-field statistics within each blob, and as 
noted above, the ordering of the blobs does not influence 
the prediction for the average extension. However, note 
that the scaling of the standard deviation, Eq. (17) differs 
from Eq. (8). We see that ctr does not depend upon 
either dw, -Dh, or w. 

For the special case of square channels, regime Ha 
has been studied under the name ‘extended de Gennes 
regime’ [12, 14, 21, 22]. Fig. 1 shows that the limits of 
this regime are more restrictive for rectangular than for 
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square channels: even if a square channel with side length 
either Du or is in regime Ha, the rectangular channel 
with side lengths Z?h and Dw might not be. 

Finally consider the ordering of length scales 

Zch < (-K < Zcw < kc <C L, (18) 

corresponding to D H <C < D w <C (Hh 4/®)^ 2 - 
This regime is denoted as lib in Fig. 1. The steps 
needed to derive the scaling relations for the extension 
are identical to those summarised above, and again lead 
to Eqs. (16,17), albeit with different prefactors than for 
regime Ha. The scaling for the extension in this regime 
was previously derived by Odijk [11]. 

Comparison with results of computer simulations. Our 
results give a qualitative explanation for the surprising 
observation of Gupta et al. [8], who found in simulations 
that while the average extension of a confined polymer 
was relatively insensitive to the aspect ratio of the con¬ 
fining channel, the variance increased with aspect ratio 
(see Fig. 5 in Ref. [8]). They performed simulations with 
a fixed channel height Z?h = 100 nm and channel widths 
in the range D w = 100 — 1000 nm, and compared these 
against simulations of square channels with matching ef¬ 
fective channel size. For their simulations, they used a 
polymer with £k = 137.4 nm and w = 18.7 nm. Thus, 
while the square channels are all approximately in regime 
Ha, the channels with fixed height cross over into regime 
lb as the aspect ratio increases significantly above unity 
(Fig. 1). While our prediction for mean extension shows 
the same scaling in these regimes, the variance is indepen- 
dent of Du and in regime Ha but increases as D^f in 
regime lb, qualitatively explaining the results of Gupta 
et al.. That the agreement is only qualitative is not sur¬ 
prising, considering that the conditions for regime lb are 
only marginally satisfied, and that the contour length of 
A-DNA which they replicate is not quite long enough to 
enter the asymptotic regime where r oc L. 

Accuracy of mean-field theory. The results that we de¬ 
rived for regimes I and II are based on Flory’s mean-field 
theory, which is thought to be correct in two dimensions 
but is known to be only approximate in three dimensions 
[33]. Using the scaling R oc L 0 - 588 in three dimensions 
would lead to modified scaling predictions for regime 
la (but not for the other regimes): R oc D^ 03 ' D^^ 3 , 
a 2 cx D-°- 37 D 2 ^ 3 , F c cx D- 1 - 70 . 

For square channels in regime Ha, we have recently 
shown both that the scalings of mean-field theory are ex¬ 
act in this regime, and that there are rigorous bounds for 
the prefactors [22]. These results were derived by map¬ 


ping the statistics of the extended de Gennes regime to a 
one-dimensional model known as the weakly self-avoiding 
random walk. Since the same mapping can be performed 
for rectangular channels, these exact results are valid 
throughout regime Ha. The scalings of Eqs. (16,17) can 
thus be rigorously proven. The rigorous bounds for the 
prefactors are included in the table in the SM [29]. As in 
regime Ha, it is in principle possible to map the statistics 
of regime lib to a solved one-dimensional model. Com¬ 
puting the exact parameters of the mapping would re¬ 
quire performing an integral over the monomer density 
profile in the Odijk regime, yet even without it the exis¬ 
tence of the mapping shows that also in this regime, the 
scalings of the extension statistics are exact. 

Free energies. Apart from the statistics of the exten¬ 
sion, the free energy of confinement F c is of experimental 
relevance, as it determines the force that must be applied 
to introduce a polymer into a channel. This free energy 
can be estimated by kT times the number of collisions 
with the walls, or 

F c /(kT) « L/l ch + L/l cw « L/l ch . (19) 

For regime la, we obtain a scaling prediction for F c 
by inserting Eq. (4) into Eq. (19), yielding F c « 
L(Ik.w/D^) 1 / 3 . For all other regimes, the free en¬ 
ergy of confinement is the same as for ideal polymers, 
F c ss L£ K /D^ in regimes lb and Ila [34], and F c « 
L(f'K-Djj) -1 / 3 in regimes Ic and lib [24]. More exact 
predictions for F c are included in the table in the SM 
[ 29 ], 

Conclusions. Recent experiments on DNA in rectangu¬ 
lar nano-channels are performed at high aspect ratio, yet 
most analytical and simulation results pertain to square 
channels. These analytical results are often applied to 
rectangular channels by matching the cross-sectional area 
of the rectangular channel to that of the square one. Our 
theory shows that this matching allows one to predict the 
average DNA extension under fairly general assumptions. 
We also show that this procedure fails to correctly predict 
the scalings of other observables, such as the variance. 
Our theory explains recent numerical results for the ex¬ 
tension variance in rectangular channels [8], and shows 
that square channels are most useful for biological appli¬ 
cations where it is beneficial that the extension variance 
is small. We expect that the results summarised here can 
be generalised to other observables, such as the statistics 
of circular DNA, knot formation, and DNA dynamics. 
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We summarise our predictions for the scalings for the 
statistics of the extension, and the free energy of confine¬ 
ment in Table I. Where numerical prefactors are known, 
they have been included in the table. In this table, we 
also include previous results for the scalings, at very 
strong confinement (regime III of Fig. 1), and for a poly¬ 
mer confined to a slit (the limit D\y —> oo). These results 
are briefly discussed below. 

1. Odijk regime 

If both Dp ip and D\y <C ip, then it is impossible 
for the polymer to turn around completely. Instead, the 
polymer must stay almost parallel to the channel axis, 
undulating slightly from side to side. In this regime the 
statistics of the extension and the free energy of confine¬ 
ment have been studied in detail both theoretically and 
numerically [SI, S2]. 

2. Backfolded Odijk regime 

Odijk [S3] has predicted the existence of a regime in¬ 
termediate between the Odijk regime and the extended 
de Gennes regime, for very slender semi-flexible poly¬ 
mers. For square channels, this regime was studied by 
Muralidhar et al. [S4], who gave it the name “backfolded 
Odijk regime”. In this regime the size of the channel is 
such that backfolds are possible but rare. Odijk defines 
the global persistence length g as the orientational cor¬ 
relation length of the corresponding ideal polymer. The 
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